Can we predict the failure point of a loaded composite material ? 
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Abstract 

As a model of composite material, the fiber bundle model has been chosen -where a bundle of fibers is subjected 
to external load and fibers have distributed thresholds. For different loading conditions, such a system shows few 
precursors which indicate that the complete failure is imminent. When external load is increased quasi-statically - 
bursts (number of failing fibers) of different sizes are produced. The burst statistics shows a robust crossover behavior 
near the failure point, around which the average burst size seems to diverge. If the load is increased by discrete steps, 
susceptibility and relaxation time diverge as failure point is approached. When the bundle is overloaded (external load 
is more than critical load) the rate of breaking shows a minimum at half way to the collapse point. The pattern and 
statistics of energy emission bursts show characteristic difference for below-critical and over-critical load levels. 
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1. Introduction 

Prediction of the failure point is a major challenge in 
various scenarios of fracture and breakdown fjl 0, lit] 
-ranging from fracturing in nano-materials to onset of 
earthquakes. A material body can tolerate certain level 
of load or force on it and beyond that level it collapses. 
If the load is increased continuously - when does the 
collapse point come ? Is there any precursor which sig- 
nals that the complete failure is imminent ? 

Fiber bundle model (FBM) has been proved flH,!!, 
0] to capture the essentials of failures in composite ma- 
terials. FBM contains a large number of fibers with sta- 
tistically distributed strength thresholds. It has simple 
geometry and clear-cut rules for how stress caused by a 
failed fibers is redistributed on intact fibers. Most im- 
portantly, this model can be solved analytically to an 
extent (For reviews, see |7[]) that is not possible for any 
other model of material-failure. The statistical proper- 
ties of FBM is well studied id §, H, and the failure 



dynamics at a constant load has been formulated 111 111 
through recursion relations which in turn explore the 
phase transition and associated critical behavior in this 
model. 

In this article we discuss how we can predict the fail- 
ure point of a loaded FBM from the available precur- 
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Figure 1 : The fiber bundle model of composite materials. 
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sors. The term precursor usually means 1 13] some prior 
indications of an upcoming incident and in current con- 
text such incident is the complete failure (collapse) of 
a fiber bundle under external load. We consider equal- 
load-sharing model, in which the load previously car- 
ried by a failed fiber is shared equally by all the remain- 
ing intact fibers [7]. The bundle consisting of N elastic 
fibers, clamped at both ends (Fig. [TJ. All the fibers obey 
Hooke's law with force constant set to unity for simplic- 
ity. Each fiber i is associated with a breakdown thresh- 
old x, for its elongation. When the length exceeds x, the 
fiber breaks immediately, and does not contribute to the 
strength of the bundle thereafter. The individual thresh- 
olds Xj are assumed to be independent random variables 
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with the same cumulative distribution function P(x) and 
a corresponding density function p(x): 

Prob(x, <x) = P(x) = f p(y) dy. (1) 
Jo 

We analyse three different loading cases: quasi-static 
loading, load increment by equal steps and overloaded 
situation. For prediction purpose it is important that 
precursors can be seen in a single system. Therefore, 
throughout this article we will present and discuss re- 
sults that can be seen in a single bundle containing 
large number of fibers. For simplicity, we consider the 
uniform distribution of fiber thresholds: P(x) = x for 
< x < 1. 
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Figure 2: The burst distributions: all bursts (squares) and bursts within 
an interval Q.9x c and x c (circles). The figure is based on a single 
bundle containing N = 10 7 fibers with uniformly distributed fiber 
thresholds within and 1 . 



2. Quasi-static loading 

The quasi-static loading is a strain controlled method, 
where at each step the whole bundle is stretched till the 
weakest fiber (among the intact ones) fails. At an elon- 
gation x per surviving fiber the total force on the bundle 
is x times the number of intact fibers. The expected or 
average force at this stage is therefore 

F(x) = Nx{l -P(x)). (2) 

The maximum F c of F(x) corresponds to the value x c 
for which dF/dx vanishes. Thus 

1 - P(x c ) - x cP (x c ) = 0; (3) 

where x c is the critical elongation value above which the 
bundle collapses. 

2.1. Burst or avalanche of failing fibers 

When a fiber fails, the stress on the intact fibers in- 
creases. This may in turn trigger further fiber failures, 
which can produce bursts (avalanches) that either lead 
to a stable situation or to breakdown of the whole bun- 
dle. A burst is usually defined as the amount or number 
(A) of simultaneous fiber failure during loading. It was 
shown in JH that the average number of burst events of 
size A follows a power law of the form 

D(A)/N = CA^ (4) 

in the limit N — > oo. Here, £ = | is the universal 
burst exponent and C is a constant. The value of £ is, 
under mild assumptions, independent of the threshold 
distribution P(x): the probability density needs to have 
a quadratic maximum somewhere within the range of 
threshold values. 



2.2. The crossover behavior 

When all the bursts are recorded for the entire failure 
process, the burst distribution follows the asymptotic 
power law D(A) oc A~ 5 ^ 2 . If we just sample bursts that 
occur near the breakdown point, a different power law is 
seen -which can be explained analytically I9i 1 1 (11 . Fig. 
|2] compares the complete burst distribution with what 
we get when we sample merely bursts in the threshold 
interval (0.9x c , x c ). 

This observation may be of practical importance, as 
it gives a criterion for the imminence of complete fail- 
ure (9[]. The bursts or avalanches can be recorded from 
outside -without disturbing the ongoing failure process. 
Therefore, any signature in burst statistics that can warn 
of imminent system failure would be very useful in the 
sense of wide scope of applicability. It is enticing to 
note the recent observation 01411 of a crossover behav- 
ior in the magnitude distribution of earthquakes before 
the largest earthquake appears. A similar crossover be- 
havior is also seen 1121 1 in the burst distribution and 
energy distribution of the fuse model which is a stan- 
dard model for studying fracture and breakdown phe- 
nomena in disordered systems. Most important is that 
this crossover signal does not hinge on observing rare 
events and is seen in a single system (see Fig. [2j. There- 
fore, such crossover signature has a strong potential to 
be used as a detection tool. 

2.3. Variation of average burst size 

We have seen that if external load is increased quasi- 
statically on a bundle of large number of fibers, bursts of 
different sizes occur during the whole breaking process 
till complete failure. One can ask - what is the average 
burst size at a particular elongation (x) value? The aver- 
age burst size is indeed a very relevant quantity that can 
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Figure 3: Average burst size A flv vs. elongation x for the same fiber 
bundle as in Fig. [2] 



ON 

S ■> 0.6 
0.4 



0.2 0.3 

X 



Figure 4: Inverse of average burst size is plotted against x for the same 
data set as in Fig. [3] A straight line can be fitted near x c from which 
one can predict the failure point. 



be measured easily during the failure process. It seems 
(Fig. [3} that average burst size( A av ) around some elon- 
gation value x goes as 



A av (x) - (x c - x) 



0.9 



(5) 



This means if we plot A^i vs. x, we should get a 
straight line which touches the X axis at x = x c . Even in 
a single system we can see this signature (Fig. [4j. 



3. Load increment by equal steps 

In the force-controlled case, load can be increased on 
a bundle by equal amount at each loading step. Here, 
the failure dynamics of the bundle can be represented 
by a recursion relation: Let N t be the number of fibers 
that survive after step f, where t indicates the number 
of stress redistribution steps. We call cr = F/N, the 
applied stress and U, = N t /N, the surviving fraction of 
total fibers. Then the effective stress after t step becomes 



u, 



and after t + 1 steps the surviving fraction of 



total fibers will be U, + \ — 1 - P(x t ). Therefore we can 
construct the following recursion relation 111 111 



U, 



t+i 



l-P(cr/U,y,U = 1. 



(6) 



At equilibrium U,+\ = U t = U*. For uniform fiber 
strength distribution, the cumulative distribution be- 
comes P(o-/U t ) = crjU,. Therefore the recursion re- 
lation becomes 



U, + i = 1 . 

U, 



(7) 



At the fixed point the above relation takes a quadratic 
form U* - U* + cr = 0, with the solutions 



t/V) = i ± (<Tc ~ Cr) 



1/2 



1 

;cr c = -. 



(8) 



One can define the breakdown susceptibility x, as the 
change of U*(<f) due to an infinitesimal increment of 
the applied stress cr: 



X = 



dU*(cr) 



dcr 



19) 



To study the failure dynamics around the critical point 
(cr — > cr e ), the recursion relation (Eq. [7} can be replaced 
by a differential equation 



dU 
~~di 



U 2 -U + cr 

u ' 



(10) 



Close to the fixed point one gets 111 111 
U,(cr) - U*(cr) * exp(-f/r), 

where the relaxation time t = | [^(cc - cry 
Therefore, near the critical point r diverges: 



TK(cr c -cry H ;6= -. 



(11) 



1/2 



(12) 



Since the susceptibility (x) and the relaxation time 
(t) follow power laws (exponent = - 1 /2) with external 
stress and both diverge at the critical stress, therefore, 
if we plot^- -2 and r~ 2 with external stress, we expect a 
linear fit near critical point and the straight lines should 
touch X axis at the critical stress. We indeed found sim- 
ilar behavior (Fig. |5]l in simulation experiments for a 
single sample. 

4. The over-loaded situation 

What happens if the initial applied load F = Ncr is 
larger than the critical load of the bundle ? The step- 
wise failure process continues and the bundle collapses 
at some step tf. If we consider the uniform distribution 
of fiber thresholds, and assume that the load is slightly 
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Figure 5: Variation of x 1 and r~ 2 with applied stress for the same 
bundle as in Fig[2] The dotted straight lines are the best linear fits near 
the critical point. 



above the critical value: cr - cr c + e, with e > 0, we can 
rewrite the recursion relation (Eq. [7]) as 



= 1 -(cr c + e)^-. 

<->t 



(13) 



The solution 111711 of the above recursion shows that 
there is a relation between the minimum of the breaking 
rate R{t) = dU t ldt (treating t as continuous variable) 
and the final step tf. 



tf = 2t . 



(14) 



At the minimum of the breaking rate, the bundle is just 
halfway to its complete collapse. Simulations on a sin- 
gle sample show that the breaking rate has a minimum 
at some value fo(e), and that for varying e the minima 
all occur at a value close to i when plotted as function 
of the scaled variable t/tf (Fig.|6]l. 



5. Energy emission bursts of loaded FBM 

It is well known that during fracturing process in 
composite materials, energy releases in the form of 
acoustic emissions and most of the cases these 

acoustic bursts follow power laws. Very recently the 
statistics of energy emission bursts has been studied 




Figure 6: Breaking rate R(t) vs. step t (upper plot) and vs. the rescaled 
step variable t/tf (lower plot) for the same fiber bundle as in Figf2] 
Different symbols are used for different excess stress levels cr - cr c : 
0.001 (circles), 0.003 (triangles) and 0.005 (squares). 



{l(h in FBM -both analytically and through numerical 
simulations. As the fibers obey Hooks law up to the 
failure, when a fiber fails at an elongation x, elastic en- 
ergy of amount jKx 2 will be released, where K is the 
force constant. Therefore for a burst of size A the corre- 
sponding energy release En can be calculated as 



A+m/n 



En 



= \ K I # 



(15) 



where x,- is the strength of failing fiber. Now at a con- 
stant applied load, if we record En at each step of load 
redistribution, it shows different pattern (Fig. [7} de- 
pending on the stress level -at critical, over-critical or 
below-critical. If we record such energy emission bursts 
separately for below-critical and over-critical levels, the 
corresponding distributions exhibit convincing differ- 
ence (Figs. 1819) . For stresses below-critical level, there 
are many many small energy bursts which are absent in 
the later case. The exponent of the energy burst distribu- 
tions are different for these two situations: -1 and —1.5 
respectively. 
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Figure 7: Energy emission En vs. step number in the same bundle as 
in Fig|2] Different symbols indicate different stress levels: at critical 
stress <t c (circles); stress cr c + 0.001 (triangles) and stress o~ c - 0.001 
(squares). 
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Figure 8: Distribution of energy emissions En in the same fiber bundle 
as in Fig|2]for stresses below critical value: cr c - 0.001 (circles), cr c - 
0.002 (triangles). The straight line has a slope —1.0, 



6. Summary and discussions 

In summary, the fiber bundle model of composite ma- 
terial shows some precursors which can help to pre- 
dict the failure point of the system under different load- 
ing situations. When external load is increased quasi- 
statically, near the failure point the exponent of burst 
distribution shows a crossover from -5/2 to -3/2 and 
the average burst size diverges. While the crossover 
signature can only warn of an imminent failure, it is 
possible to predict critical elongation value (x c ) in ad- 
vance by measuring the average burst size. If the load 
is increased by equal steps, susceptibility and relaxation 
time seem to diverge at failure stress (<x c ) following ro- 
bust power laws - from which one can predict the failure 
stress value without approaching too close to the failure 
point. When the bundle is overloaded the rate of break- 
ing has a minimum at half way to the collapse point and 
the distributions of energy emission bursts follow differ- 
ent power laws for below-critical and over-critical load 
levels. A very recent study 1I20I1 shows that for over- 
loaded case, the energy bursts also attain a minimum 
value around half way to complete collapse. 
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Figure 9: Distribution of energy emissions En in the same bundle as 
in Fig. ff]for stresses above the critical value: cr c + 0.001 (circles), 
(T c + 0.002 (triangles). The straight line has a slope —1,5. 
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